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A  Survey  of  Techniques  for  Determining 
Short-Duration  Precipitation  Rate  Statistics 


I  INTRODUCTION 

In  the  past  10  to  15  years,  there  has  been  increasing  demand  for  information 
on  short-duration  precipitation  rates  (intensities  averaged  over  intervals  of  5  min 
or  less).  Primary  interest  is  in  1-min  intensities,  commonly  referred  to  as 
instantaneous"  rates.  Knowledge  of  the  frequency-distribution  of  short-duration 
rates  is  important  to  the  design  and  operation  of  many  types  of  equipment.  Pre¬ 
cipitation,  especially  at  heavier  intensities,  attenuates  microwave  signals  of 
■nr  i'orct  systems  used  in  satellite  detection  and  tracking,  communications,  air 
traffic  control,  reconnaissance,  and  weaponry.  Erosion  due  to  rain  is  important 
to  'he  design  and  operation  of  heiieoptor  rotor  blades,  leading  edges  of  aircraft 
and  missiles,  anti  fuses  on  airborne  ordnance.  Intense  rainfall  can  cause  jet 
engines  to  malfunction  and  can  penetrate  protective  coverings  on  exposed  elec- 
•  ■■'!.  I  :  i '  .  hnnica .  materif  I. 

K'Mufaii  statistics  have  been  collected  at  thousands  of  locations  worldwide, 

‘■.r  more  than  l'i"  wars  in  many  instances.  However,  data  collection  was 
<)i  icons  toward  agricultural  anti  hydrological  purposes,  for  which  monthly,  daily, 
,es--  ■  on.  moiii  v,  :i-  and  <  -hourly  totals  were  collected.  Hates  for  '!  hr  down 
'  o  .o  in  nvsuisbit  for  many  stations  in  the  United  States,  but  for  very  few 

(is  i  iv*  o  i  *:  n.ibil  a  tint.  1*  N  iveiubor  1  !'i:2  * 
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stations  elsewhere  in  the  world.  Clock-hour  data  (totals  on  the  hour  every  hour) 
are  also  available  for  numerous  stations  in  the  United  States,  some  stations  in 
Europe,  but  only  a  few  stations  elsewhere. 

Most  rain  gauges  in  every  day  use  do  not  have  adequate  resolution  to  measure 
short-duration  rates.  Much  of  the  meager  amount  of  data  that  are  available  were 
collected  during  special  field  programs  conducted  for  brief  periods  of  time  (for 
example,  1  to  3  years).  Therefore,  researchers  have  devised  models  for  esti¬ 
mating  the  frequency  distribution  of  instantaneous  precipitation  rates  by  relating 
the  available  data  to  climatological  data  (for  example,  mean  precipitation,  mean 
temperature,  etc.  )  that  are  available  for  most  observation  sites  in  the  world. 

This  report  presents  a  review  of  techniques  for  estimating  the  frequency 
distribution  of  instantaneous  rates  at  a  point,  and  an  evaluation  of  the  limits  of 
their  applicability.  Some  authors  of  these  techniques  use  the  term  rainfall  (or 
rainfall  rate)  instead  of  precipitation.  The  use  of  either  term  in  this  report 
implies  the  inclusion  of  the  melted  equivalent  of  frozen  precipitation. 


2.  HAC.KGKOL'.NO 


There  are  three  major,  readily  available  sources  of  short-duration  precipi¬ 
tation  data.  The  largest  is  published  by  The  National  Oceanic  and  Atmospheric 
Administration  (NOAA1).  It  contains  the  highest  annual  intensities  for  intervals 

9 

of  5  to  180  min,  for  30  year's,  at  approximately  2a<>  stations  in  the  U.S.  Lin“  and 
I.er  developed  methods  to  use  these  data  for  estimating  the  distribution  of  short- 
duration  precipitation  rates  at  the  cities  for  which  they  are  available.  The  other- 
two  sources  * published  (under  contract  with  USA  I  ')  an  extensive  amount  of  data 
and  analyses  for  rates  of  4  min  or  less  at  four  networks  and  at  numerous  point 
locations.  One-minute  distributions  were  provided  for  13  of  these  locations, 


1.  National  Oceanic  and  Atmospheric  Administration  (1080)  Climatological  Data 

National  Summary,  in  annual  issues  since  lf'nl),  U.  S.  Dept.,  of  Commerce. 

2.  I. in,  S.  II.  ( 1  f > 7 7 )  Nationwide  long-term  rain  rate  statistics  and  empirical 

calculation  of  11-GHz  microwave  rain  attenuation,  Beil  System  Tech.  .J. 
5b(Nc>.  !■):  1181-  1014. 

3.  Lee,  U'.O.  V.  (IP"'1)  An  approximate  method  for  obtaining  rain  rate  statistics 

for'  use-  in  signal  attenuation  estimating,  IEEE  Trans.  Antennas  Propag. 
AP-27 (No.  3):4i>7-4  13. 

4.  Sims,  A.  L.  ,  and  Jones,  D.  M.  A.  ( 1  f 1 7 3 )  Climatology  of  Inslantam  oi.s 

Precipitation  Hates,  A  EC  K  I. -  I’ll -73 -u  17  ] ,  Pinal  Bcport,  iuntract 
E 1  f>f)28-72 -C-oil:V2  ,  Illinois  State  Water  Survey,  Urbana,  Ili.,  A 1)  7<ui~ 8-">. 

7.  Jones,  D.  M.A.,  and  Sims,  A.I..  ( 1  f *7  1 )  Climatology  of  Instantaneous 
Precipitation  Hates,  A  PC  K  I  -Til  -72 -0430,  Pinal  Repoi  I ,  <  ontract 
P  1  f 1  f !  U  8  -  C.  f  i  -  (  -OUT  ■  i ,  Illinois  State  Water  Survey,  Urbana,  Ill.,  AD  7  4  ‘ 1 1 7 ;  . 
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seven  in  the  United  States,  three  in  Europe,  and  three  in  the  Tropics.  A  poten¬ 
tially  extensive  source  of  data  for  the  U.  S.  is  discussed  by  Bodtmann  and 
Ruthroff.  They  present  a  method  for  analyzing  1-min  rainfall  rate  distributions 
from  weighing  rain  gauge  recordings,  a  large  number  of  which  are  available  from 
the  National  Weather  Service.  Results  are  presented  for  a  5 -year  period  at  2  0 
U.S.  cities. 

7 

Watson  et  al  conducted  a  survey  of  2  5  European  countries  to  determine  the 
availability  of  short-duration  rainfall  data.  Their  summary  of  the  results  reveals 
eight  countries  that  have  instantaneous  data  for  at  least  one  site,  and  19  countries 
that  have  collected  either  5-min  or  hourly  rates.  These  sources  hold  promise  for 
future  studies  of  short-duration  rates  in  Europe.  Segal  analyzed  original  rain¬ 
fall  recorder  charts  for  47  Canadian  stations.  The  results  are  in  the  form  of 
curves  of  long-term  probability  of  exceeding  instantaneous  rainfall  rates  of 
10  to  300  mm/hr  for  each  location. 

Some  researchers  (for  example,  Freeny  and  Gabbe®)  present  data  and  analysis 
of  short-duration  rates  for  individual  rain-gauge  networks.  These  provide  inter¬ 
esting  statistics  on  the  temporal  and  spatial  variability  of  short-duration  rainfall. 
Other  researchers  utilized  radar  to  measure  precipitation  rates,  but  these  are  not 
accurate  enough  for  most  applications.  The  magnitude  and  diversity  of  the  studies 
on  short-duration  precipitation  preclude  discussing  the  vast  majority  of  them  here. 
However,  Huff9 10  and  Court11  present  an  excellent  overview  of  literature  on  pre¬ 
cipitation  research.  The  emphasis  of  our  report  is  on  empirical  models  or  tech¬ 
niques  for  estimating  the  frequency  distribution  short-duration  precipitation  at  a 
point  based  on  rain-gauge  observations  over  longer  time  intervals. 


6.  Bodtmann,  W.F.,  and  Ruthroff,  C.  L.  (1970)  The  measurement  of  1  min  rain 

rates  from  weighing  raingage  recordings,  J.  Appl.  Meteorol.  15: 1 100- l  loo. 

7.  Watson,  P.A.,  Sathiaseelan,  V.,  and  Potter,  13.  (1981)  Development  of  a 

climatic  map  of  rainfall  attenuation  for  Europe,  Interim  Report  for 
European  Space  Agency,  ESTEC  Contract  No.  4  102 /7P/NL/ DG(SC'), 

Report  300,  134  pp. 

8.  Segal,  B.  (1979)  High  Intensity  Rainfall  Statistics  for  Canada,  Communications 

Research  (  enter  Report  No.  1329-E,  Dept,  of  Communications,  Ottawa, 
Canada,  124  pp. 

9.  Freeny,  A.E.,  and  Gabbe,  .1.1).  (1989)  A  statistical  description  of  intense 

rainfall,  Bell  System  Teeh.  J.  48:1789-1851. 

18.  Huff,  F.A.  (1974)  Statistics  of  precipitation,  J.  Rech,  Atmos.  ,  pp.  73-88. 

11.  Court,  A.  (197°)  Precipitation  research,  197.4-1978,  Rev.  Geophys,  Space 
Phys.  ,  Hydrology  17(No.  8):  1 165- 1 175. 
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3.  ESTIMATING  SHORT  DURATION  RATES  FROM 
LONGER  DURATION  MEASUREMENTS 


The  paucity  of  measurements  of  short-duration  rainfall  has  prompted  studies 

of  the  distribution  of  short-duration  intensities  within  a  larger  interval  of  time, 

12 

usually  1  hr.  A  very  early  study  by  Bussey  showed  how  the  instantaneous  rate 
distributed  itself  around  the  mean  hourly  rate.  He  found  that  for  Washington, 

D.  C'.  about  2  0  percent  of  an  hour  the  rate  was  a  "trace"  or  less,  35  percent  of 
the  time  the  mean  hourly  rate  was  exceeded,  and  "to  exceed  it  by  five  or  six  times 
for  a  few  minutes  was  a  fairly  common  occurrence.  "  Briggs  and  Harker  exam¬ 
ined  data  from  special  gauges  at  W  mohcombe,  England  to  obtain  the  distribution 

of  2 -min  intensities  associated  with  various  ranges  of  clock-hour  totals.  Their 

12 

results  are  in  general  agreement  with  Bussey,  and  suggest  that  they  may  be 
generally  applicable  to  showery  precipitation.  On  a  broader  scale,  Davis  and 
McMorrow^  used  actual  1-  and  4-min  rainfall  measurements  at  13  locations  to 
develop  relationships  between  clock-hour  and  short-duration  rates. 

Due  to  the  importance  of  instantaneous  rates  for  some  applications,  research¬ 
ers  have  studied  the  relationship  of  these  rates  to  5-min  rates.  For  example, 

1 5 

tfershfield  '  used  statistical  techniques  to  estimate  extreme  1-min  rates  from  the 
3 -min  rates  published  by  NOAA.  1  Pinkayan  and  Ketratanaborvornlr’  derived  an 
exponential  equation  relating  1-,  2-,  3-,  and  4-min  rates  to  5-min  rates. 

These  techniques  are  most  useful  for  the  locations  (mostly  in  Europe  and 
North  America)  for  which  clock-hour  or  other  similarly  comprehensive  precipi¬ 
tation  measurements  are  available.  For  deriving  frequency  distributions,  models 
that  utilize  data  in  the  form  of  monthly  or  annual  means  of  temperature  and  pre¬ 
cipitation,  number  of  days  with  precipitation,  etc.  ,  have  been  devised. 


12.  Bussey,  II.  E,  (Ib5D)  Microwave  attenuation  statistics  estimated  from  rainfall 

and  water  vapor  statistics,  Proc.  IRE  38:781-78.3. 

13.  Briggs,  J.  ,  and  Barker,  I.  A.  (1  pf; 2)  Estimates  of  the  duration  of  short-period 

rainfall  rates  based  on  clock-hour  values,  Meteorol.  Mag.  28:245-2 52. 

14.  Davis,  A.R.,  and  McMorrow,  D.  J.  (1275)  Stochastic  Models  for  Deriving 

Instantaneous  Precipitation  Kate  Distributions,  U.  5.  Air  Force,  Air 
Weather  Service  Tech  Report  A WS-TR -75-253,  AD  A035423,  USAFFTAC, 
Scott  A  FB,  Ill. 

15.  Hershli  Id,  I).  M.  (lr,72)  Estimating  the  extreme -value  1 -minute  rainfall, 

•  E  Appl.  Meteorol.  1 1 :2.8(i-24i'. 

1-  .  Pinka'  ui,  S.  ,  and  Ketratnnaborvorn,  T.  (ln75)  Measurement  and  analysis  of 

verv  short  duration  rainfall.  Hydrological  Sciences  -  Bulletin  -  des  Sciences 
Hydrologiques,  Bangkok,  Thailand,  XX:37-2l. 


4.  MODELS  FOR  ESTIMATING  INSTANTANEOUS 
PRECIPITATION  RATE  DISTRIBUTIONS 


Five  models  are  presented  in  this  section,  of  which  four  are  used  to  estimate 
annual  instantaneous  precipitation-rate  distributions.  The  fifth  model  estimates 
monthly  instantaneous  rate  distributions.  Three  models  require  input  of  climatic 
data,  and  two  require  only  identification  of  the  applicable  region. 

4.1  Ricc-Ilolmbcrg  Model  ( R-ll ) 

The  model  proposed  by  Rice  and  Holmberg'-'  estimates  the  percentage  of  an 
average  year  that  t-min  surface  rainfall  rates  exceed  a  given  value  (for  t  =  inter¬ 
vals  of  1  min  to  1  day).  For  1-min  rates,  the  model  uses  two  parameters  for  a 
specific  location,  the  mean  annual  precipitation  (M),  and  the  ratio  of  annual 
thunderstorm  rainfall  to  the  total  annual  rainfall  ((i).  For  t  >  1  min,  an  additional 
parameter,  the  mean  annual  number  of  days  with  precipitation  (D)  >  U.25  mm 
(o .  01  in.  )  is  also  used.  M  is  available  for  almost  all  meteorological  observation 
sites,  and  the  authors  provide  it  as  contours  on  a  world  map.  The  ratio  li  is  not 
readily  available  and  must  be  calculated  from  meteorological  data.  The  authors 
have  done  this  without  describing  the  procedure,  and  present  the  results  as  con¬ 
tours  or.  a  world  map.  The  number  of  days  with  precipitation  is  generally  avail- 
■li'U  ,  but  minimum  threshold  values  (0.  01  in,  1mm,  etc.  )  vary  from  country  to 

countrv. 

The  model  was  developed  using  data  collected  in  the  U.S.  ,  including  extreme 
short  -  duration  precipitation  published  by  NOAA  ‘  for  the  years  l«jl-lWil,  which 
w  •  extrapolated  to  estimate  1-min  rates.  W  orld  contour  maps  of  M  and  d  re¬ 
quire  interpolation  that  is  very  difficult  in  areas  where  isolines  are  closely  spaced, 
t’hese  are  major  shortcomings  of  the  model. 

1.2  Diitlnn  l)i>ii"lii'rl\  Mode!  (D-D) 

This  model  (Dutton  et  ai1^)  is  an  adaptation  of  the  R-H  model  to  Europe  and 
it  utilizes  the  same  input  parameters:  M,  ti,  and  D.  The  D-D  modification  to  the 
I! -II  mo.n  l  facilitates  error  analysis  of  the  predicted  distribution,  and  simplifies 
pro<  educes  for  estimating  ram  rates  for  a  specific  percentage  of  a  year.  The 

17.  Rice,  P.  E.  ,  and  Hoimberg,  N.R.  (11)73)  Cumulative  time  statistics  of 

surface  -  point  rainfall  rates,  IEEE  Trans.  Commun.  COM-2  l(No.  10): 

1  Id  I  -  1  IDG . 

18.  Dutton,  E.J.,  Dougherty,  H.  T.  ,  and  Martin,  R.F.,  Jr.  (1974)  Prediction  of 

European  Rainfall  and  Link  Performance  Coefficients  at  8  to  30  GHz, 
Institute  for  Telecommunication  Sciences,  U.S.  Dept,  of  Commerce, 
Boulder,  Colo.,  AD/A-000804. 


authors  also  present  a  method  for  calculating  0  for  a  specific  location  and  for 
estimating  year-to-year  variations.  Differences  in  the  distributions  between  the 
R-H  and  D-D  models  are  slight. 

Rainfall  data  for  249  stations  were  grouped  into  ten  climatic  rainfall  zones 

covering  Europe.  Values  of  M,  /3,  and  D  were  determined  for  each  station  and 

averaged  for  the  stations  in  each  zone.  These  were  used  to  construct  nomograms 

of  rainfall  rate  versus  the  percent  of  time  the  rate  is  exceeded  in  an  average  year, 

for  rate -averaging  times  of  1,  5,  30,  60,  360,  and  1440  min  (1  day),  for  each  of 

the  ten  zones.  The  distribution  for  zone  2,  which  includes  the  British  Isles  and 

western  France,  is  shown  in  Figure  1.  Contour  maps  of  M,  P,  and  D  for  Europe, 

provided  by  the  authors,  can  be  used  to  estimate  distributions  at  specific  sites. 

Contour  maps  for  Europe  of  1-min  rainfall  rates  for  1,  0.  1,  and  0.  01  percent  of 

a  year  from  the  report  are  shown  in  Figure  2. 

1  ° 

Subsequently,  Dutton  and  Dougherty  calculated  input  parameters  to  the  D-D 
model  for  305  U.  S.  weather  stations  based  on  30  years  of  data.  The  results  are 
presented  as  contour  maps  of  1-min  rainfall  rates  predicted  for  1,  0.  1,  and 
0.  01  percent  of  a  year  in  the  U.  S.  (Figure  3) 

4.3  Crane  Model  (Cr) 

2  0 

Crane  processed  seven  years  of  excessive  precipitation  data  for  15  stations 
in  New  England  and  eastern  New  York  from  the  annual  reports  published  by 
NOAA.  1  He  found  that  the  mean  and  variance  of  the  annual  extremes  were  similar 
to  the  mean  and  variance  for  24  years  at  one  of  the  stations  (Boston).  The  mean 
of  the  annual  extreme  rates  was  100  mm/hr  with  a  standard  deviation  (SD) 

35  mm/hr  in  the  15-station  sample.  The  mean  and  SD  for  Boston  were  91  mm/hr 
and  44  mm/hr.  He  felt  that  his  results  warranted  the  use  of  climatically  uniform 
regions,  and  further  developed  a  global  rain-rate  climate  model  (Crane  ).  In  it 
the  world  is  divided  into  eight  regions  (Figure  4)  based  on  total  rain  accumulation, 
and  the  number  of  thunderstorm  days  from  maps  published  by  I.andsburg.  '  The 
U.  S.  is  covered  by  five  regions;  but  one  of  them  is  further  divided  into  three  sub- 
regions  (Figure  5).  Crane  obtained  additional  guidance  from  the  Koppen  world 

19.  Dutton,  E.J.  ,  and  Dougherty,  H.T.  (1979)  Year-to-year  variability  of  rain¬ 
fall  for  microwave  applications  in  the  U.  S.A.,  IEEE  Trans,  t'ommun. 

C 'OM -27 (No.  5):82 9-837 ■ 

2n.  Crane,  R.K.  (1977)  Prediction  of  the  effects  of  rain  on  satellite  communica¬ 
tion  systems,  Proc.  IEEE  65:456-474. 

2  1.  Crane,  R.K.  (lr,8n)  Prediction  of  attenuation  bv  rain,  IEEE  Trans,  t'ommun. 
COM -2 8 (No.  9):  17  17  -  1733.  . 

22.  I.andsburg,  H.  E.  ,  Ed.  (l(l74)  World  Survey  of  Climatology,  Vols.  1-15, 
Elsivier  Pub.  (  o.  ,  Amsterdam,  The  Netherlands. 
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Figure  3.  Estimates  of  the  1-min  Rainfall  Rates  (inm/hr)  in  the  U.  S. 
Exceeded  in  a  Year:  (a)  1  Percent  of  the  Time,  (b)  0.  1  Percent  of  the 
Time,  and  £c)  For  n.  01  Percent  of  the  Time  (From  Dutton  and 
Dougherty *■  )  (Continued) 

climate  classification.  Boundaries  were  adjusted  to  accommodate  variations  in 

terrain,  predominant  storm  type  and  motion,  general  atmospheric  circulation, 

and  latitude.  Satellite  and  precipitation  frequency  data  were  used  to  extend  the 

2  3 

rain  climate  regions  over  the  oceans.  Feldman  presents  an  elaborate  expansion 
of  the  model  for  ocean  areas. 

The  measured  instantaneous  rain-rate  distributions  that  were  available  for 
each  of  the  seven  legions  and  three  sub-regions  were  pooled  to  construct  the  rain- 
rate  distributions  (Table  l).  No  data  were  available  for  Region  A,  and  the  author 
does  not  explain  how  its  distribution  was  derived. 


23.  Feldman,  N.  F.  (If?!')  Rain  Attenuation  Over  Earth-Satellite  Paths,  Final 
Report,  Contract  Nnii.i3f.-7 i'-C-Oiafi  for  Naval  Electronic  Systems 
Command,  Science  App]  i  cations,  Inc.,  A  D  A07.")3?'0. 
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Figure  t.  Global  Main  Hate  Climate  Hegions  (From  Crane^) 


Table  1.  Point  Rain  Rates  (mm/hr)  Versus  Percent  of  Year  Rain  Rate  is 
Exceeded  for  the  Regions  Shown  in  Figures  4  and  5  (From  Crane^l) 
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i.  I  Jones  Sims  Model  (J-S) 
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dunes  and  Sims  utili/a  d  1-  and  4 -min  rainfall  data  from  If'  locations  in  the 
Northern  Hemisphere  to  describe-  the'  rainfall  rate* -frequency  relationships  for 
four  of  the  five  rain  climate  regions,  worldwide,  used  by  Herbstritt.  “  This 
i.rorn'  dassification  system  (Figure  <)  does  not  include  polar  regions  for  which 
rainfall  rates  are  presumably  not  of  importance  for  most  applications.  Figure  7 
show  s  the  curves  for  the  tour  regions  derived  bv  averaging  the  frequencies  for  the 
uulr  muni  stations  w  ithiti  each  of  the  regions . 


■  Fines,  ii.  M .  A .  ,  and  ."im,  * .  I ..  (i'Mti)  *  luiiatniogy  of  mstanunn  ous  ramlal! 

rates,  d  .  App  1 .  ' 1 1  U  c •  i .  I  V  (No .  n  I :  i  I  d  ■  -  1  1  4  ■  >. 

!  1 1  :  i  -s t;a  tt ,  !< .  !  . .  i  !  ■  -7  1  I  c.  i -  i  1  *  i  ■  i i a  ti  t  mis ;  r  ra t ions  lor  I  ioim  ;  i  ■  >a  ti  1 1  i ' m-  - 

ontltK  n*al  I'.-'..  Id  ; m  !  No.  i.-7  ."I,  !•  •  -'era.  1  -  an  :m:r.  n,- hate  '  mu in  i>  - 

-ion,  ihii-'i  ol  i.ia  1  iic  :  ini.uif.  i,  A  i-duimton,  I ).  t  .  ,  •  pi  . 


Figure  fi.  World  Rain  Climates  Used  by  Herbstritt' 


The  obvious  shortcomings  of  the  J-S  model  are  the  mixing  of  1-  and  4 -min 
data,  ana  the  use  of  extensive  regions.  Resolution  is  poor  when  information  tor 
individual  locations  is  required. 


4.5  Lenhard  Model  (Le) 

The  only  method  we  found  that  can  be  used  to  estimate  monthly  rainfall  rate 

->  p 

statistics  was  developed  by  Lenhard.  ~  He  used  data  on  1-min  rainfall  rates  for 

5  4 

13  stations,  published  by  Jones  and  Sims  and  Sims  and  Jones,  and  developed  a 
model  for  estimating  instantaneous  rainfall  rates  equalled  or  exceeded  2.  f),  1.0, 
0.  5,  0.1,  n.  05,  and  0.  01  percent  of  the  time  during  a  month  from  rei  lily  avail¬ 
able  climatological  data.  The  model  is  of  the  form 


A  +  B  I  +  C  T 
P  P  P 


(1) 


where  Fl^  is  the  precipitation  rate  (mm, /min)  equalled  or  exceeded  p  percent  of 
the  time  during  the  month,  I  is  a  precipitation  index  (the  ratio  of  mean  monthly 
precipitation  and  the  number  of  rainy  days),  and  T(°F)  is  the  monthly  mean  tem¬ 
perature.  A  ,  b  ,  and  C  are  coefficients  that  depend  upon  the  exceedance  prob¬ 
ability  ((>).  bach  coefficient  is  given  by 


K  =  <•  +  r  In  p  f  y  !n“  p 


(2) 


where  K  is  the  coefficient  (A  ,  or  C  )  and  <>,  ;J,  and  7  are  least-squares 
, ■•caress :or.  coeffi e icnts. 

Separate  sets  of  eciuations  were  derived  for  pi  ecipitation  indices  base-:  or. 
days  w  itn  n.  nj  in.  ,  1.  1  m.  .  or  1  mm  or  mo:  <  to  define  a  rainy  daw  Another 

threshold  value  to  define  a  rainy  day  is  a  trace,  but  different  os  between  the 
nmr.be"  of  days  w  ith  this  amount  and  1.  ■)]  in.  or  more  w ere  found  to  be  slight, 
fable  gives  tile  •  oef: i .-icnts  to  be  used  ir.  J.q.  (ft >  lor  each  of  tilt  three  dofimtiun.- 
01  a  -lav  with  ram.  A  Iso,  eot  ffieienti;  arc  present  n.  separately  as  derived  using 
■ia'a  !  ah  1  •  stations  in  the  study  .  iron-  me  ten  «  xtra -tropical  "’ations,  and 
Iron  if-  'in.  *  trap:  a!  stations.  \s  an  examine.  to  estimate  the  raml'a.L  rau 
e\  •  1  •  .  .  |  ••••<  n‘  •>:  *r,<  ti me  :: !  !  ’  ■  ns  ii.  .* u i  .  ;.  (  1 '  b.  ■  01m  s 


q  .in'.-  :. 


For  an  index  of  5.8  mm  per  rain  day  based  on  days  with  1  mm  or  more,  and  a 
mean  temperature  of  67.5°F,  the  equation  yields  an  estimate  of  0.86  mm/min. 


Table  ‘>.  Coefficients  to  Produce  Estimating  Equations  for  the  Lenhard  Model 
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I.cnhard  cautions  that  the  model  cannot  be  extended  below  the  input  data, 

(that  is,  mean  temperature  <22°I  and  index  values  <2  mm/day).  In  fact,  some 
internal  inconsistencies  were  found  at  index  values  below  10  mm/day  and  tem¬ 
peratures  below  4()°F.  Rates  should  decrease  with  an  increase  in  the  probability 
of  >  xrccdanrc.  To  determine  if  the  equations  are  applicable  to  a  location  that  is 
i  ithcr  very  cold  or  very  arid  or  both  is  to  calculate  rates  for  all  exceedance 
probabilities.  If  they  are  consistent  with  each  other,  they  are  acceptable.  If 
they  are  consistent  but  a  negative  rate  shows  up  at  the  higher  exceedance  prob¬ 
abilities  (for  example,  1  or  2  percent),  it  could  indicate  that  it  doesn't  rain  that 
often  in  that  month. 

l.fi  Comparisons 

I’lie  D-I>  model  is  a  modification  of  the  R-H  model,  and  actual  differences  in 
the  distributions  calculated  using  them  a;  e  slight  (Dutton  et  al^).  One  of  the 
input  parameters  for  both  models  is  the  ratio  of  thunderstorm  rain  to  total  rain 
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at  a  location.  Dutton  et  al  presents  an  expression  to  calculate  this  ratio,  but 

it  requires  the  highest  monthly  precipitation  observed  in  30  consecutive  years  of 

data.  Few  stations  outside  of  North  America  and  Europe  have  recorded  such  data. 

17 

However,  a  world  map  of  the  ratio  is  provided  by  Rice  and  Holmberg. 

The  Cr  and  J-S  models  provide  instantaneous  rate  distributions  for  a  given 
percent  of  a  year  for  regions  delineated  for  this  purpose.  Distinct  advantages 
of  the  Cr  model  are  the  greater  number  of  regions  used  and  the  greater  clarity 
of  regional  boundaries  (see  Figures  4  to  6).  A  descriptive  comparison  of  the 
models  is  provided  in  Table  3. 

A  comparison  of  estimates  using  the  R-H,  D-D,  Cr,  and  J-S  models  with 
observed  values  for  U.  S.  cities  is  given  in  Table  4.  The  observed  values  were 
taken  from  the  distributions  derived  from  five  years  of  original  records  at  each 
city  by  Bodtmann  and  Ruthroff.  These  data  were  not  used  in  the  development  of 
any  of  the  models.  R-H  model  estimates  were  made  using  M  from  long-term 
climatic  data  instead  of  the  world  contour  map  for  M  provided  with  the  model. 

U-D  model  estimates  were  made  using  their  contour  map  for  the  U.  S.  [Figure  3(c)]  . 
The  best  estimates,  as  indicated  in  Table  4  by  the  root-mean-square  departures, 
were  made  using  the  Cr  model,  followed  by  the  D-D,  R-H,  and  J-S  models. 

The  Le  model  was  not  compared  to  the  others  in  Table  4  because  it  produces 
estimates  of  instantaneous  rates  on  a  monthly  basis.  Although  these  statistics  are 
probably  more  useful  than  annual  distributions  for  many  applications,  no  other 
monthly  or  seasonal  global  models  were  found.  At  AFGL,  development  of  an 
improved  model  for  estimating  monthly  disti  ibutions  is  being  completed.  A  report 
on  the  model  will  be  released  in  the'  near  future. 


->.  1)1  Kl\l\(.  INSTANT ANKOl  S  PRECIPITATION  RATE 

STATISTICS  FROM  CLOCK  IIOl  REV  DISTR IIJLTIONS 

Four  models  are  presented  in  this  section.  Two  use  clock-hour  data  to  derive 
Mu-  distribution  of  instantaneous  rates,  and  two  can  be  used  to  derive  clock-hour 
•  list  ributions. 


•7.1  M.mI.  F  tor  I >cri\ inu:  Instantaneous  I  )i -I  dilutions 

Davis  and  M.-Morrov  ^  derived  tables  of  -'look-hour  rale  intervals  versus 
instantaneous  rate  intervals  for  sb.  location-:,  and  versus  4 -in  in  rates  for  seven 
other  locations,  worldwide.  Hie  tables  give  the  percent  contribution  of  the  short- 
duration  rate  w  ithin  <  ach  of  nine  spe  i  died  intervals  for  each  i  lock-hour  rate 
interval  Tin  location  for  which  the  short -duration  precipitation  rate  distribution 
is  rcquiri  d  must  be  compared  v  ith  tile  stations  in  Mu  report.  Tin-  tables  l’ox-  the 
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Table  3.  Comparison  of  Models  for  Estimating  Instantaneous  Precipitation  Rates 
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Table  4.  Comparison  of  Model  Estimates  of  Instantaneous  Rates  With  Those  Ob¬ 
served  at  Nine  U.  S.  Cities.  Departures  from  observed  values  are  given  in 
parenthesis 


Rate  for  0.  01  Percent  of  an 

Avg  Year  (mm/hr) 

Location 

Observed 

R-H 

D-D 

Cr 

J-S 

Atlanta,  Ga. 

69 

93  (+2  4) 

92  (+2  3) 

63(-6) 

100(+3 1) 

Bismarck,  N.D. 

41 

30(-il) 

35(-6) 

37 (-4) 

6  5  (+2  4) 

Buffalo,  N.Y. 

48 

63  (+5) 

48(0) 

43*  (-6) 

63(1-17) 

Columbus,  Ohio 

66 

62  (+6 ) 

5  5  ( - 1  > 

4  9 (—7 ) 

6  5  (+9) 

Dallas,  Tex. 

73 

85(-l2) 

7  5  (+2 ) 

8 6’-  (-17) 

100(+27) 

Memphis,  Trim. 

64 

93  (+2  9) 

80(+16) 

G  3  ( —  1 ) 

100(436) 

Miami,  Ela. 

111 

115(4-4) 

11  >(+4) 

f>8(- 13) 

K'O(-ll) 

Milwaukee,  Wis. 

:>■> 

of)  (+4) 

a.)  (+3) 

43  (-9) 

68(4-13) 

Newark,  N.J. 

)(l 

43 (-2 ) 

t:n(+  j  m) 

4:,(- 1 ) 

6  3(4-15) 

Root-mean -squat 

e  departures 
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Estimates  of  instantaneous  rates  from  observed  clock-hour  rates  using  the 
Davis  and  McMorrow  tables,  and  from  Huschke's  table  that  combines  all  the 
Davis  and  McMorrow  data,  were  determined  for  0.  01  percent  of  a  year  at  the 
nine  U.  S.  cities  in  Table  4.  These  estimates,  provided  by  the  USAF  Environ¬ 
mental  Technical  Application  Center  (ETAC),  are  given  in  Table  5  along  with 
observed  values  for  each  city  from  Table  4.  ETAC  made  the  Davis  and  McMorrow 
estimates  using  the  Urbana,  Ill. ,  table  for  all  locations  in  Table  5,  except 
Memphis  and  Miaini,  for  which  the  Franklin,  N.  C. ,  and  Miami  tables,  respec¬ 
tively,  were  used.  Tables  to  estimate  1-min  rates  are  also  available  for  Majuro 
Atoll,  Marshall  Islands;  Woody  Island,  Alaska;  and  Island  Beach,  N.  J.  The  Island 

Beach  table  would  seem  the  logical  choice  to  make  the  estimate  for  Newark  in 
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Table  5.  However,  Davis  and  McMorrow  incorporate  an  internal  check  to 
determine  the  applicability  of  the  analog  model  to  the  rain-rate  statistics  for  the 
station  being  modeled.  They  indicate  that  the  mean  precipitation  at  a  location 
correlated  best  with  the  lower  35th  percentile  integrated  rainfall  amount  calculated 
from  the  estimated  instantaneous  rate  distribution.  The  35th  percentile  value  for 
Newark  is  110.5  cm  (43.5  in.),  and  118.4  cm  (46.  6  in. )  using  the  Urbana  and 
Island  Reach  models,  respectively.  Since  the  actual  annual  average  precipitation 


Table  5.  Comparison  of  Instantaneous  Bates  Estimated  From  Clock -hour 
Hates  With  Those  Observed  at  Nine  U.  S.  Cities.  Departures  from  ob¬ 
served  values  are  given  in  parentheses 


Location 

Rate  for  0.  01  Percent  of  an  Avg  Year  (mm/hr) 

Observed 

Davis  ft  McMorrow 

Huschke 

Atlanta,  Ga. 

69 

7l(+2) 

7  2  (+3) 

Bismarck,  N.  D. 

41 

38(-3) 

3  5  ( —  6 ) 

Buffalo,  N.Y. 

48 

47  (- 1 ) 

44  (-4 ) 

Columbus,  Ohio 

56 

54 (-2 ) 

52 (-4) 

Dallas,  Tex. 

73 

6P(-4) 

72  (-1) 

Memphis,  Tenn. 

64 

65(+l) 

7  9  (+ 1 5 ) 

Miami,  Fla. 

111 

10P(-2) 

1 02  ( -  9 ) 

Milwaukee,  Wis. 

52 

5  4  (+2) 

5 1  ( - 1 ) 

Newark,  N.J. 

50 

5  8  (+8) 

54  (+4) 

Hoot-mean-square  departures 

3.  4 

6.  7 

for  Newark  is  105.4  cm  (41.  5  in.  )  the  Urbana  model  was  used.  The  rate  estimated 
f<-  •  0.01  percent  of  the  year  using  the  Island  Beach  model  is  60  mm/hr. 

The  root-mean-square  departures  from  the  observed  values  using  the  Davis 
and  McMorrow,  and  Huschke  models,  given  in  Table  5,  are  smaller  than  for  the 
models  in  Table  4.  These  results  justify  the  use  of  clock-hour  data,  when  avail¬ 
able,  to  estimate  instantaneous  rate  distributions.  When  not  available,  the  models 
in  the  next  section  can  be  used  to  estimate  clock-hour  distributions. 

5.2  Models  for  Deriving  Clock-Hour  Distributions 

2  8 

Using  clock-hour  data  for  four  U.S.  sites,  Russak  and  Easley  found  good 
linear  correlations  between  a  climatological  index  (Cl)  and  the  number  of  hours 
per  year  that  specified  hourly  rates  are  exceeded.  Cl  is  simply  the  ratio  of  the 
mean  annual  precipitation,  in  inches,  to  the  mean  annual  number  of  rainy  days. 

They  determined  that  the  regression  lines  of  Cl  versus  the  hours  per  year  (11^) 
the  rate  (r)  is  exceeded  comprised  a  family  of  lines  that  could  be  expressed  by 
the  relationship 

H  =  m(CI  -  a)  ,  (") 

where  m  is  the  slope  of  the  regression  line  for  rate  r,  and  a  is  the  x-axis  inter¬ 
cept  of  the  regression  Line  for  rate  r.  It  should  be  noted  that  five  threshold  values 
to  define  a  rainy  day  are  in  general  international  use;  trace,  0.  1  mm,  1  mm, 
n.  01  in.  ,  and  0.  1  in. 

Russak  and  Easley  provide  curves  of  m  anil  a  versus  precipitation  rate,  so  the 
user  needs  only  the  mean  annual  precipitation  and  number  of  rainy  days  at  a  loca¬ 
tion  to  use  the  model.  Bv  determining  the  number  of  hours  per  vear  of  rainfall  in 
the  interval  between  two  rates  for  a  sufficient  number  of  intervals,  it  is  possible 
to  construct,  a  cumulative  frequency  distribution  for  a  location.  This  was  done  bv 
the  authors  for  Naples,  Italy,  and  the  results  compare  favorably  with  the  frequency 
distribution  from  ac  tual  clock-hourly  data.  They  note,  however,  that  the  method 
is  not  valid  in  certain  climatic  regions.  For  example,  erroneous  results  are  likely 
in  arid  locations,  areas  with  strong  monsoonal  or  orographic  control,  or  regions 
where  convective  precipitation  is  very  predominant.  The  model  accuracy  probably 
suffers  because  it  does  not  address  differences  in  the  threshold  rainfall  amount 
used  to  define  a  rainy  day. 


28.  Russak,  S.  L.  ,  and  Easley,  I.W.  (1P58)  A  practical  method  for  estimating 

rainfall  rate  frequencies  directly  from  climatic  data.  Bull.  Am.  Meteorol. 
Son.  35(No.  0):4G0-472. 


Another  model  for  estimating  the  annual  distribution  of  clock -hour  rates  at  a 
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location  was  developed  by  Winner.  He  used  clock-hour  data  from  123  U. S. 
stations  to  develop  nomograms  that  can  be  used  to  determine  the  number  of  hours 
per  year  that  equal  or  exceed  hourly  rates  of  trace,  0.01,  0.  02,  0.  10,  0.25,  0.  50, 
1  and  2  in.  /hr.  One  of  the  input  parameters  is  the  ratio  of  the  mean  annual  pre¬ 
cipitation  to  the  average  annual  number  of  days  with  precipitation  5:0.  01  in. 
Another  parameter  in  the  model  is  called  a  "moisture  index"  for  which  information 
on  potential  evapotranspiration  is  needed.  The  author  suggests  two  sources  for 
this  data. 

This  model  was  intended  to  overcome  the  geographical  limitations  of  the 
Russak  and  Easley  model.  However,  results  compared  poorly  with  independent 
data  from  the  tropics.  Therefore,  the  author  decided  to  develop  a  different 
method  for  the  tropics.  Using  data  from  32  stations  in  the  Panama  Canal  zone, 
he  obtained  best  results  from  a  linear  correlation  of  the  annual  rainfall  with  the 
number  of  clock-hours  at  the  various  rates.  The  resulting  relationships  for 
estimating  the  mean  annual  number  of  clock-hours  with  precipitation  equal  to  or 
greater  than  0.25,  0.50,  0.75,  and  1.00  in.  for  the  tropics  are  provided  in 
Figure  8. 

A  deficiency  in  Winner's  approach  is  the  threshold  value  of  0.  01  in.  for  a 
rainy  day.  Although  this  value  is  used  in  the  U.  S.  ,  it  is  not  common  in  most 
other  countries.  Another  shortcoming  is  the  use  of  only  U.  S.  and  Panama  data. 


ft.  OTIII  K  CONSIDERATIONS 

Accuracy  of  instruments  used  to  record  rainfall  is  affected  by  the  height  and 

size  of  the  collector,  the  wind  speed,  and  freezing  conditions.  Gauges  also  have 
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mechanical  limitations,  especially  at  high  rates.  Grayman  and  Eagleson  report 
that  these  result  in  5  to  10  percent  error  in  total  rainfall  catch.  Larger  errors 
are  likely  in  measurements  of  instantaneous  rates  derived  from  standard  record¬ 
ing  tipping  or  weighing  bu.  ket  gages  This  is  due  to  poor  resolution  of  the  record¬ 
ing.--  down  to  ]  min.  Special  observation  programs  have  generally  utilized  high¬ 
speed  r<  i  orders  to  minimize  errors  in  determining  short-duration  rates;  however, 
•In -si  data  arc  limited  in  quantity  and  generally  available  for  only  short  periods 
(that  is,  1  to  years). 

-  ■  A  mm  r,  I).  <  .  (l',r:;)  <  Tiinatologieal  Estimates  of  C’loek-Hour  Rainfall  Rates. 
U,  S.  Air  Loi  re,  Air  U  father  Service,  Technical  Report  202. 

( ;•  -ivii.aii,  A.M.,  and  Eagleson,  P.  ?.  (inVo)  A  Review  of  the  Accuracy  of 
Rada  r  and  Raitigagt  s  tor  Precipitation  Measurement.,  Hydrodynamics 
E.*l  ora  mi' v  K<  port  No.  11  I  >ept.  of  Civil  Engineering,  MIT,  Cambridge, 
Mass. 
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Year-to-year  variations  in  the  distribution  of  instantaneous  rates  observed 

at  a  location  can  be  quite  large.  The  models  presented  in  this  report  represent 

the  distribution  to  be  expected  on  average  at  a  point.  These  may  not  be  adequate 

for  many  design  problems  for  which  a  low  risk  is  important.  Annual  variations 
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are  discussed  for  the  U.  S.  by  Dutton  and  Dougherty,  for  Canada  by  Segal,  and 
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for  Europe  by  Dutton  et  al,  and  Dougherty  and  Dutton. 

7.  CONCLUSIONS 

When  instaritaneous  precipitation-rate  statistics  are  needed  for  a  specific 

point  location,  it  is  best  to  derive  such  distributions  from  actual  measurements 

made  continuously  over  a  period  of  "many"  years.  Section  2  of  this  report  refers 

to  sources  of  such  distributions  for  many  locations  in  North  America,  (for  exam- 

pie,  Segal,  and  Bodtmann  and  Ruthroff  )  and  potential  sources  of  data  for  the 

U.  S.  (NOAA,  *  and  Bodtmann  and  Ruthroff*’)  and  Europe  (Watson1).  Distributions 

f”om  observations  taken  during  special  studies  for  short  periods,  1  to  3  years, 

5  *  1 

for  13  stations  worldwide  were  derived  by  Jones  and  Sims  and  Sims  and  Jones. 

For  locations  where  eioek-hour  data  are  available,  good  estimates  of  instantaneous 
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rates  can  be  made  using  models  by  Davis  and  McMorrow  or  Huschko.  "  Data 
for  short  periods,  that  is,  less  than  five  years,  should  be  used  with  caution  since 
they  might  not  give  a  true  indication  of  the  distribution  of  rates  during  an  average 
year. 

Where  distributions  from  actual  observations  are  not  available,  they  can  he 
determined  from  the  five  models  discussed  in  Section  4.  They  require  the  input 
of  climatic  data  that  arc  available  for  many  observation  sites  worldwide.  The 
appropriate  model  to  use  depends,  in  part,  on  the  application  and,  in  part,  on  the 
climatic  data  available  for  model  input.  Hie  model  developed  by  l.c  nhard"  is  the 
only  one  that  can  be  eised  to  estim.b.e  instantaneous  rates  versus  frequency  on  a 
monthly  basis  for  a  specific  location.  The  other  four  models  can  !><  used  to  deter¬ 
mine  instantaneous  rates  versus  percent  of  a  year.  Two  of  their,  pre  sent  distribu¬ 
tions  that  represent  all  the'  individual  locations  within  regions.  (  'f  the  se  ,  the 
model  developed  by  Crane'''  is  recommended  lie-cause  it  has  a  ere  ate  i  mmiie  i  oi 
regions  with  more  < -learly-clc line  d  boundaries  than  the  othe  r  mei.ie  1  :e  ve  Vepe-  I  by 

•>4  ... 

Jones  and  buns.-  I  wo  nmdcls  can  be'  use-el  tee  estimate  tile  distributin'!  .  >  l  instan¬ 
taneous  rates  for  a  year  based  on  cJimatu  data  at  individual  ineatmns,  howe  ver 
one-  is  actually  a  refine-ment  eiver  the  oilier.  Hie  original  nioiU  i,  Pi  ve  iupi  >i  by 

3  1.  Dougherty,  H.  T.  ,  ami  Dutton,  I.J.  ( 1 : 1 T :  )  Estimating  >e  ar-lo-ye  availabil¬ 
ity  of  rainfall  for  microwave  applies tions,  II  I  11  Frans,  e  etn.m an. 

C  OM  -2fi(No.  8 ) :  132  1-1324. 


Klee  and  Hoitnberg.  was  modified  by  Dutton  et  at10'  1  and  applied  to  a  large 
number  of  locations  in  Europe  and  the  U.S.  to  produce  maps  ot  rate  versus  per¬ 
cent  of  a  year.  Two  models  for  estimating  clock-hour  distributions  are  also 
presented,  but  the  use  of  these  estimates  with  either  the  Davis  and  McMorrow, 
or  Huschke  models  are  not  likeiv  to  produce  better  results  than  the  Crane,  Dutton- 
Dougherty,  or  ltice-Holmberg  models. 

The  precision  of  distributions  of  instantaneous  rates,  whether  derived  from 
actual  data  or  from  an  empirical  model,  is  limited  by  instrument  error  and  poor 
resolution  of  recordings  down  to  l  min.  The  combination  of  errors  is  likely  in 
excess  of  10  percent.  Furthermore,  the  models  we  have  presented  estimate 
distributions  to  be  expected  on  average  at  a  point,  but  year-to-year  variations  in 
instantaneous  distributions  at  a  location  can  be  quite  large.  This  is  discussed  in 
Section  (■>  and  should  be  considered  in  design  problems  for  which  a  low  risk  is 
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